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ABSTRACT 

New and previously published observations of the bright eclipsing binary VV Ori- 
onis are analyzed. We present new radial velocities and interstellar reddening measure- 
ments from high-resolution spectra of this detached, short-period (P=1.48 d) binary. 
We discuss the validity of prior claims for the existence of a third body and show that 
our new velocities and light curve solution cast doubt on them. The components of 
VV Ori are shown to be a Bl V primary with a mass Mi = 10.9±0.1Mq and a radius 
Ri = 4.98 ± O.O2i? and a B4.5 V secondary with a mass M 2 = 4.09 ± O.O5M and 
a radius R 2 = 2.41 ± O.OTi? . 
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1 INTRODUCTION 

VV Orionis (HD 36695) is a bright (V=5.4), double-lined 
eclipsing binary consisting of main sequence B-type stars 
in a detached configuration. The inclination is high enough 
■ to produce complete eclipses, enabling the very accurate 
determination of fundamental parameters of the system. 
Since the discovery of its photometric variability over a 
century ago, the system has been the subject of sev- 
eral spectroscopic and photometri c studies as discussed by 

ISarma fc Vivekananda Raol (|l995). Unfortunately, the de- 
rived absolute parameters in these previous studies vary 
considerably and this led us to obtain additional spectro- 
scopic and photometric data on the system as par t of ou r 
program on early-t ype binaries (viz. iTerrell. et alJ (2003), 

iTerrell. et all (120051) ). 

Datin g back to th e first comprehensive spectroscopic 
stu dy bvlDaniell l lSld) , spec t roscopic investigators , includ - 
ing IStruve fc Luvtenl (119491) IBeltrami fc GaleottH J1969D . 

lDuerbeckTi^75l) and lPoppen lm)93T) ""nave concluded that a 
third star orbits the eclipsing pair. The evidence for claims of 
this third body rests solely on higher than expected resid- 
uals in fits to the radial velocities of the primary as the 
third body's lines are not seen in the system's spectrum. 

iDaniell lll915l) offered the third body h ypothesis "with great 
reserve" and IStruve fc Luvtenl II1949T) pointed out that the 
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observed variation in the systemic vel ocity did not nec essar- 
ily imply that the system was a triple. iDuerbeckl dl975l) ana- 
lyzed his radial velocitie s alo ng with those of IDaniell 
IStruve fc Luvtenl Jl949l) and IBeltrami fc Galeottil (Il969l) to 
arrive at orbital parameters for a third body orbiting with 
a period of about 119 days. A goal of the current study was 
to investigate this third body hypothesis. 



2 OBSERVATIONS 
2.1 Photometry 

VV Ori was observed with a 0.25m Schmidt-Cassegrain tele- 
scope, Santa Barbara Instrument Group ST-7XE CCD cam- 
era and standard BVIc filters. The main goal of the new 
photometric observations was to cover at least one primary 
eclipse for the purpose of extending the baseline of eclipse 
timings so that any small changes in the period could be 
detected. Given the extensive photometric datasets previ- 
ously published, we did not attempt to observe complete 
light curves for the system. 

The CCD images were reduced using IRAF. Bias and 
dark frame subtraction and flatfielding were done in the 
usual way. PyRAF scripts automated the frame calibra- 
tions and aperture phometry which was performed using 
the IRAF DAOFIND and PHOT routines. The comparison 
star for the differential photometry was HD 36779. 
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Figure 1. The appearance of the He I 5876 A line in the spectrum 
of VV Ori. The line components originating from the primary and 
secondary star are labeled " p" and " s" , respectively, and 6 is the 
orbital phase. The interstellar sodium doublet and the telluric 
H2O absorption lines longward of 5885 A are also marked. 

2.2 Spectroscopy 

High resolution, high signal-to-noise spectra of VV Ori were 
obtained with the Echelle+CCD spectrograph of the 1.82m 
telescope operated in Asiago by the INAF Astronomical Ob- 
servatory of Padova. The CCD is an E2V thinned, back 
illuminated lKxlK chip of 13.5 /J.m pixel size. The slit 
was always aligned with the parallactic angle and its pro- 
jected width was 1.5 arcsec. The spectra cover the range 
from atmospheric cut-off around 3250A to 7100A at a re- 
solving power 30,000. The accuracy of the zero point of the 
wavelength scale is derived on each observed spectrum by 
cross-correlation against a synthetic template of the rich tel- 
luric absorption spectrum at wavelengths 6860-7050A, 6270- 
6310A, 5880-5950A and found to be affected by no system- 
atic shift larger than 0.2 km sec -1 , consistent with the high 
efficiency of the inst rument in deriving accura te radial ve- 
locities evaluated bv lMunari fc Lattanzil jl992l) . The expo- 
sure times for the spectra were 5 minutes (0.002 in orbital 
phase), eliminating any concerns about phase smearing ef- 
fects on the radial velocities. In all, 31 spectra were obtained 
over the period from November, 2004 to February, 2006. 

Figure shows the appearance of the He 1 5876A line 
near quadratures. The weakness of the secondary's lines il- 
lustrates the need for high-resolution, high signal-to-noise 
observations when trying to measure the radial velocity of 
the secondary. Because of the weakness of the secondary's 
lines, its radial velocities could only be accurately measured 
near quadratures when the blending with the primary's lines 
was minimized. 



2.2.1 Reddening 

The Nal lines have two interstellar components with ra- 
dial velocities of +9.3T0.2 and +26.1T0.1 km sec -1 . The 
FWHM of both components is 14.4 km/sec, which is domi- 
nated by the instrumental PSF. The equivalent widths of the 
two components are 0.066±0.003 and 0.160±0.006 A for the 
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Figure 2. The interstellar line of Nal D2 with the two-component 
Gaussian fit used in the reddening determination. The thick line 
is the observed profile and the thin line shows the two-component 
fit. The dashed lines show the individual components of the fit. 



Nal D2 line at 5890A, and 0.040T0.002 and 0.130±0.005 A 
for th e Nai Dl line at 5896 A. Using the lMunari fc Zwitteil 
lll997t) calibration on Nal line at 5890 A this corresponds to 
a reddening of E(B-V)=0.078±0.004. 

The Tycho B T - V t = -0.20 ±0.01 corresponds, ac- 
cording to lBesselll fcOOOft . to a Johnson B - V = 0.17 ±0.01. 
We denote an intrinsic (i.e., unreddened) color by (B — V)' 
and an observed color by (B — V) with subscripts 1, 2 and C 
denoting the primary component, the secondary component 
and the composite value for the binary respectively. With 
an assumed intrinsic color for the primary, (B — V)'\, and 
the magnitude differences between the two components in 
each filter, Ams and Amy, from our light curve solution, 
we can compute the intrinsic composite color of the binary, 
(B - V)'a, as 



(B - V)' c = (B — V)'x - 2.5 log 



1 1 jQ(-Q.4xAm B ) 
1 + 10(-°- 4xA "H') 



Assum ing a value of (B - V)\ = -0.27 for the B1V pri- 
mary JStraizvslll992T) and the Am_g and Amy values of 
2.58 and 2.49 respectively from the light curve solution, we 
find (B — V)'o = —0.26. This result, combined with the ob- 
served (B - V)c, yields a reddening of E(B - V) = 0.09, 
in good agreement with the results from the analysis of 
the interstellar Nal lines. These valu es also agree well with 
the value of 0.07 measure d by both iDuerbeckl jl975h and 
IChambliss fc Leund ifi?: 



2.2.2 Radial velocities 

Radial velocity measurements for the primary have been de- 
rived from measurements of seven hydrogen Balmer lines (a, 
/3, 7, 8, e, H8, H9) and eight He I lines (7065, 5876, 5016, 
4472, 4388, 4026, 4009, 3820 A). 

Given the weakness of its spectral features (cf. Figure 
0, the measurement of the radial velocity of the secondary 
star has been possible only in a few spectra, those with the 
highest S/N near quadrature, and for only a subset of the 
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Table 1. Heliocentric radial velocities of VV Ori. 



HJD 


Primary R.V. 


Secondary R.V. 


(mid-exposure) 


(km sec ) 


(km sec ) 


2453314.5470 


-88.29 




2453373.4771 


97.01 




2453373.4804 


93.64 




2453373.4828 


95.12 




2453392.4492 


140.27 


-241.88 


2453392.4538 


142.14 


-241.59 


2453392.4581 


142.12 


-246.59 


2453393.4475 


-77.80 


288.75 


2453393.4521 


-76.09 


275.65 


2453393.4565 


-72.99 


269.44 


2453394.3750 


52.35 




2453394.3794 


47.78 




2453394.3841 


47.40 




2453411.3610 


-37.56 




2453411.3689 


-34.81 




2453411.3767 


-32.80 




2453411.4244 


-11.41 




2453411.4323 


-7.48 




2453411.4400 


-7.59 




2453412.3884 


-58.46 




2453412.3963 


-57.64 




2453412.4044 


-63.07 




2453775.2399 


-59.49 




2453777.3065 


144.07 


-274.19 


2453777.3125 


144.42 


-282.39 


2453777.3389 


142.22 


-279.10 


2453780.2507 


147.63 




2453780.2551 


147.51 




2453780.3954 


118.14 


-234.42 


2453780.4003 


114.54 


-230.26 


2453780.4050 


112.87 


-226.48 



He I lines. Table 2 gives the radial velocities of both compo- 
nents. 



3 DATA ANALYSIS 

We performed a sim ultaneous analysis of the 
IChambliss fc Leund (I1982T) UBV and vby photometry, 
our UBV photometry and our radial velocities with the 
2003 version of the Wilson-D e vinne y ijWilson fc Devinnevl 
119711: IWilsonl Il979t IWilsonl Il990l: hereafter. WD) pro- 
gram. We original l y inc luded the u observations of 
IChambliss fc Leund l|1982r) but found, as they did, that 
we could not get satisfactory fits to those data. Inspection 
of the light curves shows that the system is detached so 
we used WD mode 2 in the solution. The parameters 
adjusted in the simultaneous solution were the semi-major 
axis of the relative orbit (a), the binary centre of mass 
radial velocity (V^), orbital inclination (i), secondary mean 
effective temperature (T2), modified surface potentials 
of the components (Qi and mass ratio (g), and the 

bandpass-specific luminosity of the primary (Li). We 
used time as the independent variable and adjusted the 
orbital period (P), and the reference epoch (HJDo). No 
evidence for an orbital eccentricity was found in light curve 
fits so we fixed it at zero for the final solution. Certain 
parameters, such as the bolometric albedos and gravity 



brightening exponents, were held fixed at their expected 
theoretical values. The logarit hmic limb darkenin g law 
was used with coefficients from IVan Hammel £993) . The 
rotational velocities of the primary (167 ± 4 km sec" 1 ) and 
the secondary (83 ± 8 km sec -1 ) are consistent with the 
assumption of synchronous rotation and that was assumed 
in the solution. The mean effective temperature of the 
primary was set to 261 99 K based on the unreddened B — V 
and the calibration o f lFlowe"rl (I1996T) . Kurucz atmosphere 
models (|Kuruc3[l993) were used for both stars. Data set 
weights were determined by the scatter of the observations. 

Given the frequent claims of the existence of a third 
body in the system, we also adjusted third light but we could 
not fin d statistically signific a nt va lues for any of the light 
curves. IChambliss fc Leund <ll982T) also found no indica- 
tion o f third light in their data.lS a rma fe Vivekan anda Raj 
(1995) anal yzed the UBV data o f lDuerbeckl Il975l) and the 
H a data of IChambliss fc Davanl jl987l) and also found no 
evidence of third lighT Tciiarnbnssrifi984) . however, did find 
a third light contribution but did not give error estimates 
on which to judge the significance of the result. 

Table [5] shows the results of the si multaneous solution 
and F igures El and |H show the fits to the lChambliss fc Leund 
( 1982) data. Figure |5] shows the fits to the radial velocities. 
The primary appears to have absolute dimensions consis- 
tent with a slightly evolved main seque nce star of its spec - 
tral type. Using a stellar evolutio n code ijHa n^etal] Jl994l) , 
lEggletonl Jl973l) . lEggletonl lll972ft . lEggletonl jl97im kindly 
supplied to us by P. Eggleton, we find that the VV Ori pri- 
mary reaches its current radius at an age of 8.3 Myr, with 
the error in the radius leading to an uncertainty in the age of 
0.2 Myr. At that age, the models predict a secondary radius 
of 2.44 Rq in reasonably good agreement with our observed 
value of 2.41 g i ven t he uncertainties in the masses. 
IChambliss fc Leund (I1982T) found that the secondary was 
much smaller than expected for the mass they derived. Our 
derived mass is considerably smaller and our radius slightly 
larger than their values, easing this mass-radius discrepancy. 
Table|3]shows our results for the component masses and radii 
compared wi th pre viously published values. 

IPopperl Jl993l) estimates the spectral type of the sec- 
ondary as B5 V. Our derived effective temperatur e results 
in a B 4.5 c lassification according to the tables of IStraizvsl 
lll992ft and iBertone. et alJ <2004h . The radi us of 2.41 R ? 
is consistent with an unevolved B4 V star JStraizvsl I1992T) 
and the surface gravity is also consistent with that spectral 
type. The secondary's bolometric mag nitude is -1.57 and 
that value is in good agreement with the lStraizvsl il992f) ta- 
bles for a B4.5 V star. Experiments with composite reference 
spectra and the strength of the He I 5876 line both point to 
a B4.5 V spectral classification. Based on these factors, we 
estimate that the secondary's spectral type is B4.5 V. 

Assuming a ratio of total to selective absorp tion of 
R = 3 .1 and using the bolometric corrections from iFlowerl 
jl99_6j) , we find that the distance to VV Ori is 388 ± 30 par- 
sees with the error dominated by the uncertainty in the effec- 
tive temperature of the primary. Unfortunately, the Hippar- 
cos parallax is not very precise, yielding a la distance range 
of 390pc to 1060 pc. Our distance places the system in the 
Orion OBI as sociation cons istent with the dist ance to the 
la subgroup fela auwl Il9 64[) measured b y Ide Zeeuw. et alJ 
l|l999|) . ISarma fc Vivekananda Raol (119951) estimate that the 
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Table 3. Comparison of measurements of the masses and radii of the components of VV Ori. 



Parameter 


Current Study 


Popper (1993) 


Chambliss & Leung (1982) 


Duerbeck (1975) 


Eaton (1975) 


Ri (Rq) 


4.98 ± 0.03 


5.0 ±0.2 


4.94 


4.47 


5.0 


R 2 (Rq) 


2.41 ±0.02 


2.45 ± 0.1 


2.36 


2.21 


2.3 


Afi (M Q ) 


10.9 ±0.2 


10.8 ± 0.6 


10.2 


7.60 


9.3 


M 2 (M ) 


4.09 ± 0.08 


4.5 ±0.8 


4.5 


3.42 


4.3 



Table 2. Parameters of VV Ori 



Parameter 


Value 


a 


13.49 ± 0.05Rq 


Vy 


23.5 ± 0.3 km sec -1 


i 


85°.9±0°.2 


Ti 


26199 K 


T 2 


16073 ± 42 K 


f2i 


3.152 ±0.004 


n 2 


3.433 ±0.008 


<? 


0.376 ± 0.001 


HJD 


2452500.2065 ± 0.0003 


P 


1.48537423 ± 0.00000005 days 


L 1 /(L 1 +L 2 ) U 


0.9367 ± 0.0005 


Li/{L\ +L 2 ) B 


0.9150 ± 0.0005 


Li/(Li + L 2 ) v 


0.9083 ± 0.0005 


Li/(ii +L 2 ) Ic 


0.8999 ± 0.0014 


L 1 /(L 1 +L 2 ) V 


0.9165 ± 0.0005 


Li/(Li +L 2 ) b 


0.9111 ±0.0005 


Li/(Li + L 2 ) y 


0.9081 ± 0.0004 


Ri 


4.98 ± 0.02 R 


n 2 


2.41 ± 0.01 R 


Mi 


10.9 ±0.1 M 


M 2 


4.09 ± 0.05 M Q 


Li 


10600 ± 1600 L 


L 2 


350 ± 54 L 


log gi 


4.08 ± 0.06 


log gi 


4.29 ±0.06 



Quoted errors are formal la errors from the solution. Luminosity 
errors are estimates based on an estimated uncertainty of 1000 K 
in the effective temperature of the primary. The log g values are 
in CGS units. 

age of VV Ori is 10 ± 1 Myr. This value, along with the age 
found from our evolutionary calculations, supports the mem- 
bership of VV Ori in the la subgroup of Orion OBI whose 
age was measured to be 11.4 ± 1.9 Myr. The barycentric 
radial velocity we measure, 23.5 ±0.5 km sec -1 , agrees well 
with the average radial vel o city o f the Ori OBI members 
found bv iMorrell fc Levatol (^L991). At the galactic latitude 
of —18° and our distance, VV Ori is 110 pc from the galac- 
tic plane, within the scale height of the galactic thin disk. 
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Figure 3. The uvby photometry from Chambliss and Leung 
(1982) and the computed light curves. The u data were not used 
in the final solution and are shown only to illustrate the nature 
of the fitting difficulties. 



1.4 




0.7 -J , , 1 . 1 , 1 . 1 . , , 1 1 

-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

Phase 

Figure 4. The UBV photometry from Chambliss and Leung 
(1982) and the computed light curves. 



4 THE PUTATIVE TERTIARY COMPONENT 

iDaniell lll915t) was the first to suggest the existence of a ter- 
tiary component in VV Ori. Larger than expected residuals 
from the fit to his radial velocity observations of the pri- 
mary led him to examine the possibility of a third body and 
he gave orbital parameters for it "with great reserve ... and 
not w ith any idea that t hey a re even approximately true." 
Later. Istruve fc Luvtenl lll949l) also noted systematic trends 



in the residuals to the fits of their velocities but did not feel 
justified in solving for the orbit of a third body, stating that 
"... we cannot be certain that th e slow variation n ecessarily 
means that the system is tr iple." | Duerbeckl ( 1975) an alyzed 
the ve locit y residuals for thelDa niell lll915h . lStruve fc Luvtenl 
Jl949l) and lBeltrami fc Galeottil lll970D data sets, as well as 
his own, concluding that a third body orbit with an period 
of about 119 days fit the data best. Assuming an inclination 
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Figure 5. The radial velocities from the current study and the 
computed curves. 



of 90° for the orbit of the third body, he estimated that its 
spectral type was mid-to-late A and its luminosity was 13-14 

We repeated our light and velocity c u rve analys i s after 
add ing the velocities from [ Popper! l)l993T) . IDuerbeckl (|lflTnT l 
and IStruve fc Luvtenl Jl94Sf) to study the behaviour of the 
velocity residuals for each radial velocity data set. All radial 
velocities were given equal weight in the solution. Figure HJ 
shows the velocity residuals for the four s ets of data phase d 
with the third body ephemeris given bv IDuerbeckl Jl9T5l) . 
While the older data appear to show systematic trends, our 
more precise data do not, but we, unfortunately, do not have 
full coverage of the purported third body orbit. 

One intriguing result of the analysis is a statistically 
significant value for the time derivative of the orbital pe- 
riod of the binary. The solution converged on a value of 
P = -3.6 x 1CT 10 ± 1.3 x KT 10 , perhaps indicating ac- 
cretion of the primary's wind by the secondary. We do not 
attach much significance to this result, its reality having 
to be assessed against additional precise radial velocities 
and accurate timing of eclipses to be secured in the years 
to come. In the case of conservative mass transfer between 
the two components, the mass exchange would amount to 
2 x 10 -7 Mq yr~ x . This material would be easily ionized by 
the radiation field of the Bl V primary, and should give rise 
to emission lines easily detectable in high resolution spectra. 
Ours do not show any hint of emission lines originating in 
VV Ori 

Another way to test for t he ex istence of the third body, 
as pointed out bv IDuerbeckl il975t) . is through its effect on 
the O-C diagram for times of minimum of the binary. Unfor- 
tunately, VV Ori has not been a frequent target of observers 
of times of minimum. A concerted effort to observe as many 
times of minimum as possible over a couple of observing sea- 
sons could prove valuable in answering this critical question. 
IDuerbeckl (Il975t) predicts that the double amplitude of the 
variation would be 0.0014 days, a measurement within reach 
of a small telescope equipped with a CCD camera. 

Given the large scatter in the older radial velocity data 
sets, the lack of any evidence for third light in the light 
curves and the lack of systematic trends in the residuals of 



Figure 6. The radial velocity residuals for the simultaneous fit to 
several data sets, phased with the third body ephemeris given by 
Ducrbcck (1975). Triangles are from the observations of Struve 
and Luyten(1949), squares from Duerbeck (1975) and open circles 
from Popper (1993). Data from the current study are shown with 
filled circles. 



our radial velocities, we see no strong evidence to support 
claims of the existence of a third body in VV Ori. More ex- 
tensive radial velocity and time of minimum studies should 
be able to resolve the question with greater certainty. 



5 CONCLUSIONS 



VV Ori consists of a Bl V primary with M i = 10.9±0.1M ( 
and 

M 2 



© 



Ri = 4.98 ± 0.02%, and a B4.5 V secondary with 
= 4.09 ± 0.05M m and R 2 = 2.41 ± 0.01%. The de- 



rived distance of 388 ± 30 pc and barycentric velocity of 
23.5 ± 0.3 km sec -1 are consistent with the system's mem- 
bership in the Orion OBI association. We find no evidence 
of third light in the photometric data and the new radial 
velocities, although somewhat limited in temporal coverage, 
show no evidence of the third body discussed by previous 
authors. 
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